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We consider cosmological consequences of a heavy axino, decaying to the neutraUno in R-parity 
conserving models. The importance and influence of the axino decay on the resultant abundance of 
neutraUno dark matter depends on the lifetime and the energy density of axino. For a high reheating 
temperature after inflation, copiously produced axinos dominate the energy density of the universe 
and its decay produces a large amount of entropy. As a bonus, we obtain that the upper bound on 
the reheating temperature after inflation via gravitino decay can be moderated, because the entropy 
production by the axino decay more or less dilutes the gravitinos. 

PACS numbers: 95.35.-|-d, 98.80.Cq 



I. AXINO 

Neutralino, if it is the lightest supersymmetric particle 
(LSP) in R-parity conserving models, is a natural candi- 
date for dark matter. Because of the TeV scale sparti- 
cle interactions, the thermal history of neutralinos allows 
the neutralino dark matter possibility. But, imposing a 
solution of the strong CP problem, the thermal history 
involves contributions from the additional sector. 

The strong CP problem is naturally solved by introduc- 
ing a very light axion a. Most probably, it appears when 
the Peccei-Quinn (PQ) symmetry is broken at a scale of 
fa- Below the PQ scale, the effective axion interaction 
with gluons is 



C = 



9 s 



aFF, 



(1) 



where gs is the strong coupling constant The PQ 
scale is constrained by the astrophysical and cosmological 
considerations in the narrow window 10^° GeV < fa ^ 
IQi^GeVH. 

TeV scale supersymmetry (SUSY) suggests axino a, 
the superpartner of axion, around the electroweak scale 
in the gravity mediation scenario. Here, we consider the 
effects of heavy axinos in cosmology. The axino cosmol- 
ogy depends crucially on the axino decoupling tempera- 
ture 0], 



Ta-dcp = 10" GeV 



V/^O-iy (2) 
1012 GeV/ la'' ^ ' 



where as = gl/^ir. 
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The axion supermultiplet includes axion, saxion (the 
scalar partner) and axino. Both saxion and axino masses 
are split from the almost vanishing axion mass if SUSY is 
broken. The precise value of the axino mass depends on 
the model, specified by the SUSY breaking sector and the 
mediation sector to the axion supermultiplet Q . In prin- 
ciple, the axion supermultiplet is independent from the 
observable sector in which case we may take the axino 
mass as a free parameter of order from keV to a value 
much larger than the gravitino mass Q. Light ax- 
inos can be a dark matter (DM) candidate, which has 
been studied extensively 0, B [a]- Heavy axinos, how- 
ever, cannot be the LSP and can decay to the LSP plus 
light particles. This heavy axino decay to neutralino was 
considered in the literature @ where the neutralino relic 
density was not considered seriously. Some considered 
the axino as the next LSP decaying to the gravitino LSP 
in the gauge mediated SUSY breaking scenario [l^l . Re- 
cently, supersymmetric axion models were studied with 
an emphasis on saxion pd| . where the heavy axino pos- 
sibility was also considered briefly [l2| • 

In this paper, we present a more or less complete cos- 
mological analysis of a heavy axino with mass in the 
TeV region so that it is heavier than the LSP neutralino. 
Compared with the saxion study, the heavy axino study 
probes SUSY directly because the axino carries the odd 
R- charge. 

If kinematically allowed, a heavy axino decays pre- 
dominantly to a gluino and a gluon and the gluino sub- 
sequently decays, finally producing the neutralino LSP. 
However, if axino is lighter than gluino, this axino to 
gluino decay is forbidden and the axino predominantly 
decays to a neutralino and a photon. The neutralinos 
from axino decay can annihilate in the cosmos if the 
neutralino number density, n, is large enough so that 
n{cFannv) > H whcrc aann IS the annihilation cross sec- 
tion. In this case, the neutralino abundance is modified, 
which is obtained by solving the Boltzmann equation. 
For this to be compatible with the observed DM den- 
sity, we find the required thermally-averaged cross sec- 
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tion of neutralino, {aannv), to be around 10~®GeV~^. 
The Higgsino-hke neutrahno can give this kind of large 
cross section. 

The thermally produced 0(100 GeV) gravitinos after 
inflation induce severe problems on the light element 
abundances, which restricts the reheating temperature 
to Tr < 10^"'' GeV. We find that the entropy produc- 
tion from the heavy axino decay can dilute the primordial 
gravitinos. For axino mass smaller than gluino mass and 
for a low PQ scale, fa ~ 10^" GeV, we find as shown 
below that there is no gravitino problem for the reheat- 
ing temperature up to the axino decoupling temperature 

-^a— dcp- 

A heavy axino leads to different physical consequences 
depending on its mass being greater or smaller than the 
gluino mass uig. If axino is heavier than gluino, it de- 
cays dominantly to a gluino plus a gluon. If axino is 
lighter than gluino, it decays dominantly to a 6-ino-like 
neutralino and a photon and, if kinematically allowed, to 
a neutralino and a Z-boson. Thus, the axino decay width 
is given by, 

(i) 7715 > nig 



r{d ^ g + g) ^ 



128^3 /2 



1--4 , (3) 



(ii) rriy^ < rria < rrig: 



r(a Xt+l) 



em^aXi7 '"'a 
128^3 J2 



(4) 



with Caxi-y = {CaYY / COS 9w)Zx^B Hi where Z^.b is the 
6-ino fraction of the i-th neutralino and is a Weinberg 
mixing angle. Hereafter we will use CaVY — Z^.b ~ 1 
for simplicity. For the case (ii), as one can see, the life- 
time can be easily longer than 0.1 second. If the life- 
time is shorter than about 0.1 second, the axino decay 
does not harm the standard big bang nucleosynthesis 
(BBN). Note that for Eq. ^ the axino lifetime is about 
3.3 X 10-^s (a^/O.l)"^ (/a/10" GeV)^ (ma/1 TeV)"^ 
Right after the axino decay, the temperature of radia- 

,1/4 



tion Td = V^aMp /{n^g ' is given by 



Td = 1.4 GeV 



70 



1/4 



3 X lO-'^sec 



1/2 



(5) 



where Mp is the reduced Planck mass. For the case that 
axinos dominate the energy density of the universe be- 
fore they decay, we can regard Td as the second reheating 
temperature due to the axino decay. This temperature is 
marginal when we discuss its effect on neutralino DM as 
we discuss below. This is because, for the axino lifetime 
of 10~^ second, this temperature is comparable to the 
typical neutralino freeze-out temperature Tfr sa m^/25 
[H]. If Td > Tfr or equivalently ra < 0(10"^) sec- 
ond, the axino decay has no effect on the neutralino DM 



abundance, because the neutralino produced by the ax- 
ino decay also could reach the thermal equilibrium. From 
Eq. Q , we can see that this would be the case for a heav- 
ier axino or for a lower PQ scale. 

The axino abundance depends on the thermal history 
of the early universe after inflation. Another relevant 
temperature we introduce is the reheating temperature 
after inflation Tr. So, the temperatures we introduce are 

Ta-dc-p = axino decoupling temperature 

Tr = reheating temperature after inflation 

Tfj. = neutralino freeze — out temperature (6) 

Ta^rad = axino — radiation equality temperature 

Td = radiation temperature right after a decay. 

For Tr > Ta_dcp, axinos were in the thermal equi- 
librium and the axino number density is comparable to 
the photon number density. Below the axino decoupling 
temperature, the axino number in the comoving volume 
is conserved if there is no new physics below Ta_dcp; then 
the axino abundance is given by 



na 

^ a 

S 



135C(3) ga 
87r4 g 

* (^a— dcp) 



for 



7r, > Ta-dcp, 
(7) 



where ^(3) ~ 1.202, ga is the degrees of freedom of axino 
and s = ^^g*sT^ is the entropy density. 

For Tp{ < Ta_dcp, axinos could not be in thermal equi- 
librium after inflation. Nevertheless, they are regener- 
ated by thermal scattering and by decays of gluinos, 
squarks and neutralinos in the thermal plasma. For a 
high enough reheating temperature, T 10** GeV, the 
hard thermal loop approximation is good enough and the 
axino abundance can be approximated as [l^ 



Ya = 2.0 X 10-\g^ln 



1.108 

9s 



10" GeV y 

Ta J 



104 GeV 
(8) 



We can see that the axino abundance might be large in 
terms of Ya going up to 0(0.002), strongly depending on 
Tr and fa- Hence, if many axinos are produced after 
inflation, axinos might dominate the energy density of 
the universe before they decay. 

Let us first examine the condition for axino domina- 
tion. For T < ma, axino becomes nonrelativistic. When 
axinos are nonrelativistic, the energy densities of radia- 
tion and axinos are given by 



PR = ' P^^\T<m~^ = maYaS, 



(9) 



respectively. From Eq. ([9]), we find the axino-radiation 
equality temperature, Ta=rad, 



(10) 
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If this equality occurs before axino decays as shown as the 
inequaUty, then axino can dominate the universe. On the 
other hand, for ImaFa < 7b, axinos never dominate the 
energy density of the universe before they decay. In this 
regard, we find the lowest reheating temperature, T^"™, 
above which axinos can dominate the universe before 
they decay, by solving the equality |TOai^a(T™") = ^b- 
In case axinos dominate the universe, the entropy pro- 
duction by axino decay dilutes the previously existing 
number densities. The ratio of the entropy per comoving 
volume before and after the axino decay is [ll] 



Sq 



3Td 



for Ta=rad > Tu- Thus, from Eqs. 
entropy ratio can be 



(11) 

and (HI]), the 



So 



1.3 



VlTeV/ 



2GeV\ 



0.002 



(12) 



This is shown in Figs. [T] and [2] as magenta lines. Above 
the solid line denoted as r = 1, axinos can dominate the 
universe and can produce additional entropy. 

Similarly, in case that an axino can decay only into a 
neutralino and a photon [Case (ii)], Td looks to be able 
to become as small as of several MeV and Sf/So might 
be as large as 0(10^ - 10^). 



II. RELIC DENSITY OF NEUTRALINO 

Neutralino, which was in the thermal equilibrium in 
the early universe, decouples and freezes out when the an- 
nihilation rate becomes smaller than the Hubble param- 
eter. The freeze-out temperature Tfr is normally given 
by m^/2b, e.g. 4GeV for 100 GeV neutralino. As men- 
tioned above, the neutralino relic density is not affected 
by axino for Td > Tfr- 

Hence, we consider Td < Tfr below and the following 
arguments are valid only for this case. Another crucial 
issue in this case is whether the produced ncutralinos 
from the axino decay would annihilate again or not. If 
the neutralino number density from axino decay is too 
large, i.e. for 



90 

9 

.9* 



1/2 



1 



1 



4(cru) A/pTc 



(13) 



the produced ncutralinos would annihilate more. 



not annihilate. From Eqs. ([TO]) and ([T3)) . this occurs if 



1/2 



{(Tv) < 



AMpTuY^iTu) MIpTl 



90 



1/2 



(14) 



is satisfied. The equality in (fH)) is estimated using 
Eq. (|12[) assuming that the neutralino abundance from 
axino decay is much bigger than that from thermal freeze- 
out. This condition is not likely to be satisfied for the 
range of our ma and Td. 



This leads us to the case that Eq. (|T4)) is not satisfied, 
i.e. neutralinos produced (from axino decay and thermal 
freeze-out) annihilate. The final abundance is obtained 
solving the Boltzmann equation. 



ann ^re 



i)n 



(15) 



where aann is the annihilation cross section of two neu- 
tralinos, Vrei is their relative velocity and (. . .) is the 



thermal average. With Y^^ i 
tion is 

^Y'x / \^^2 

where s is the entropy density, and s oc 
the radiation dominated era. 



^/s, the evolution cqua- 



(16) 



T3 oc t-3/2 



Applying the sudden decay approximation, after the 
decay of axinos, the radiation and neutralino densities 
can be easily estimated. The radiation dominates the 
universe soon and the neutralinos start to annihilate. We 
can solve this evolution equation from the time after ax- 
ino decay, Td, to later times. The final abundance Y^ 
at temperature T after neutralino annihilation can be 
expressed as 



Y-\T) =Y-\Tu)-{a. 
c,Y-\Tu) + 



ann'^rel) \ ^ 
-iVrel) s{Ty)) 



H{Tu) 



(17) 



H{Tn) 



Here, Y^ (Td) is the sum of neutralino densities from ther- 
mal freeze-out and from the axino decay, right after axino 
decay but before neutralino annihilation. 



k,(Td) ^r/' + x,''^""^. 



(18) 



When axino dominates before the neutralino freeze-out, 
T'a=rad > Tfr, then the freezed-out neutralino relic den- 
sity, which is different from the standard one in the ra- 
diation dominated universe, can be obtained as 



Case of axino domination 



For Td < Tfr, let us consider the cases of axino domi- 
nation. First, let us consider the case that axinos would 
dominate the universe but the produced neutralino can- 



^fr^m7r%ry^3fgD_Y^' 
A{av)MpTf^2\gfr, 



Td 



(19) 



where qd = g*{Tu), 9fr = 9*{Tf^) and Tf^ is the freeze- 
out temperature of the neutralino in the nonrelativistic 
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FIG. 1: The Tr vs. ma plot for fa = 10^° GeV. The re- 
gion Tr > Ta-dcp is above the dashed blue line (horizon- 
tal). The axino lifetime greater than 0.1 sec is denoted by 
the red shaded region in the left side. The blue shaded region 
in the right side is where axino decays before neutralino de- 
couples (To > Tfr). The magenta lines (horizontal) are the 
contours of the entropy increase, r = Sf/So- Above r — 1 
lines axinos dominate the universe before they decay. The 
green lines (vertical) denote the (aannVrei) in units of GeV"'^ 
which are used to give the right amount of neutralino relic 
density by Eq. ((22}. We use neutralino and gluino masses as 
m 

mass, the green lines move to the right 



X — 100 GeV and nig = 2TeV. For a heavier neutralino 



axino matter dominated universe. Otherwise, i.e. for 
7a=rad < Tfr we obtain 



^ A{<7v)MpTfr 



(20) 



The relic density of neutralinos from the axino decay 
can be given by the axino abundance suppressed by the 
entropy production, 



ydecay _ yi ^ ^2. 



(21) 



where Ya is the axino abundance just before axino decay. 

When the number density of neutralinos from 
the axino's sudden decay is much larger, so that 
"■V (Jn ) {cFgnn Vrpi ) ^ H{T'£,), thcu wc obtaiu the result 
of [l5|. In this limit, the relic density can be approxi- 
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FIG. 2: The same as Fig. [T]but with fa = 10^^ GeV 
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mated and simplified as 



i0.14 



90 



1/2 



7r^5* (7b) 
10-^ GeV~^ 



100 GeV 
2GcV\ 



(22) 



where we normalized Td for an axino decay to gluino and 
gluon. The relic density of neutralino is proportional to 
the neutralino mass which is different from the standard 
one where Q,~^h? oc Xf = rriy^/Tfr ^ 25. 



B. Case of axino nondomination 

In the case that axinos never dominate the universe, 
Ta=rad < 7b, the resultant neutralino LSP abundance 
is given by Eqs. (fTT]) and p8|) . but with Y^^ given by 
eq. (Eni) and Y^'""'y = ^a, in which Tfr is the freeze out 
temperature of the neutralino in the radiation dominated 
universe, instead of Eqs. (fTO|) and (PT|) . 

If Eq. p3)) is satisfied, as in the previous case, we ob- 
tain the same final neutralino relic abundance as that 
given in Eq. (P^ . This is evident if we think that the 
final neutralino relic density after re-annihilation is de- 
termined from the Hubble parameter at the time of axino 
decay. 

On the other hand, if Eq. is not satisfied, wc sim- 



5 




mg (GeV) rrig (GeV) 

FIG. 3: The same as Fig. [T]but with = 300 GeV. FIG. 4: The same as Fig. [2] but with mg = 5 TeV. 



ply obtain 



(23) 



This is possible when the axino abundance is too small, 
which is the case for low reheating temperature Tr < 
0(10^ GeV). For a large annihilation cross section, the 
first term becomes negligible and the second term is dom- 
inant. In this case, the neutralino number density is given 
by thermally produced axino's. Hence, in this case also, 
we may have a chance to measure the reheating temper- 
ature as in the axino DM recently pointed out 
by Choi et. al. 0. 

In Figs. [T] and m we show the contour lines (green lines) 
of (av) which gives the neutralino relic density suitable 
for the DM in the universe. As can be seen in the fig- 
ures, the high annihilation cross section of neutralinos 
of order (av) > 10~* GeV is necessary. Here, we used 
the neutralino mass of 100 GeV. If we increased the neu- 
tralino mass, then the (av) increases accordingly since 
neutralino relic density is proportional to the neutralino 
mass as shown in Eq. (P^ . which is plotted in Fig. [3] 



with 



= 300 GeV. For completeness we show the 



plots with different gluino mass in Fig. |4l where we used 



mg=5 TeV. 



proportional to the reheating temperature Tr. If the 
gravitino is not the LSP then it can decay. Since the in- 
teraction of gravitino is suppressed by the Planck mass, 
the lifetime is much longer than that of axino, usually 
falling in the BBN or post-BBN era. The decay products 
of gravitino can change the abundances of light elements 
and gives the severe constraint on the reheating temper- 
ature, Tr < lOS-'^GeV for - TeV [13, 111- 

For a high reheating temperature in our scenario, the 
axino dominates the universe and excessive entropy is 
produced as shown in Eq. (fT2|) . which may soften the 
gravitino problem. This entropy production is propor- 
tional to the reheating temperature, which is of the same 
form as the gravitino production. Therefore, as Tr in- 
creases, the increased number of gravitinos is diluted by 
the increased entropy from the axino decay, which ren- 
ders the gravitino abundance independent of the reheat- 
ing temperature after axino enters into the axino domi- 
nated phase in the universe. This independence is true 
only when Tr is smaller than the axino decoupling tem- 
perature Ta-dcp- 

The abundance of gravitinos after axino decay can 
be expressed as 



S 



(24) 



III. SOFTENING THE GRAVITINO PROBLEM 

In the early universe, gravitinos are produced after re- 
heating. It is known that the gravitino abundance is 



where Y^{Ty{) is the thermally produced gravitino abun- 
dance after reheating. T^*" is the lowest reheating tem- 
perature above which the axino dominates: r = 1, de- 
noted by the solid magenta line in the figures. This is 
vahd for T^™ < Tr < Ta-dcp- 
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So for fa = lOi^GeV, this hne is 
10® GeV, where the gravitino problem is almost solved 
anyway. However, for = 10^^ GeV, T^"*" = 
around 10® GeV, the gravitino problem still exists, 
though it is tolerated in our heavy axino scenario. 



almost T™*" 



IV. CONCLUSION 

We have discussed the heavy axino possibility so that 
the heavy axino decay leads to a reasonable DM density 
of the LSP iieutralino and enough radiation to dilute the 
gravitinos. For this scenario to be realized, we need a 
large annihilation cross section of neutralinos (av) . Such 
a neutraliiio is Higgsiiio-like, in which case the cross sec- 
tion with nuclei for the direct DM search can be as large 
as (Tsi — 10"'^"®-' pb, which can be probed in future DM 
search experiments. 
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matter dominated universe, Eq. p9|) for completeness. 
Although one may find the same discussion in Ref. p^ . 
some dependence is missing there and we will correct 
it. Just after the decoupling of neutralino, from the rel- 
evant Boltzmann equation psp . we obtain 



1 



1 



(av)2 



(A.l) 



for t ^ tfr- In the matter dominated universe 
tfr = 2/iHfr, and the freeze-out temperature is Ty^ — 

— - — - Here, in the decaying 



-m^ In 



matter dominated universe, the entropy density per co- 
moving volume is not conserved. Hence, wc used a'^ in- 
stead of the entropy density in the standard calculation. 
After the reheating by the decaying matter is completed, 
when the temperature is T/, of course, we can use Y as 
usual. We find 



Y\ 



T=Tf 



3H 



a{Tf)J 2{av)s{Tjy 



fr 



(A.2) 



after the decay. Recalling g*{T)T^ cx H and a{Tf (x 
g^{T)^T^ in the decaying matter dominated universe, we 
obtain 



Y 



90 



1/2 



1 



9*iTf) \ 

g*s{Tf)J \Tr^g*{Tfr) J 4{av)MpTfr 



T. 



1/2 



2 \TfrJ \g*{Tsr), 

In text, we used g*{Tf) = g^,s(Tf) and Tf = T-q. 



(A.3) 



APPENDIX: NEUTRALINO FREEZE-OUT 
DENSITY IN DECAYING PARTICLE 
DOMINATED COSMOLOGY 

In this appendix, we derive an approximate relic den- 
sity of thermal neutralinos decoupled in the decaying 
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